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Following a solution generating technique introduced recently by one of us, we transform the 
Einstein static Universe into a two - fold infinity class of physically acceptable exact perfect fluid 
solutions of Einstein's equations. Whereas the entire class of solutions can be considered as gen- 
eralizations of the familiar Tolman IV solution, no member of the class can be written explicitly 
in isotropic coordinates. Further, except for a set of measure zero, no member of the class can be 
written explicitly in curvature coordinates either. 
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I. INTRODUCTION 

Over the last decade or so the study of exact spheri- 
cally symmetric perfect fluid solutions of Einstein's equa- 
tions has been reinvigorated [IJ. Recently, one of us [2] 
transformed the Einstein static Universe into an exact 
physically acceptable perfect fluid solution of Einstein's 
equations and examined the properties of the solution in 
familiar curvature coordinates. The emphasis there was 
on the importance of the integration constant that the 
generating technique introduced. Here we transform the 
Einstein static Universe into a two - fold infinity of so- 
lutions all of which may be considered generalizations of 
the familiar Tolman IV solution. (This solution, and its 
variants, is useful for the study of internal properties of 
neutron stars Q.) It is shown that no member of the 
class can be written explicitly in isotropic coordinates 
and, except for a set of measure zero, no member can be 
written explicitly in curvature coordinates either. 

When looking for an exact solution of Einstein's equa- 
tions, the search most often starts from a familiar set of 
coordinates. The present paper serves to emphasize the 
important point that physically interesting solutions can 
be automatically excluded by this approach. 



II. SIMPLIFIED GENERATING TECHNIQUE 

Following '^l we consider spacetimes O that can be 
written in the form Q 



dr2 



1 - 



2M(r) 



+ r^dn^^dt^\ (1) 



where dVi^ = dO^ -f sin(6')^(i(/i^ and M is constructed so 
as to make O a perfect fluid 5]. That is, 



M = 



fb{r)ei^^^''^'^''dr -f Ci 



a{r)dr 



(2) 



where = Ci is a constant, 



r(l + 2rx') 



and 



~ ^ '2.r{r(x' -x'^)-x) 
1 + 2rx 



(4) 



It is important to note that O is not a conformal trans- 
formation of the seed metric given by the case x — ^ 
(due to the restrictions on M). Moreover, M is not the 
effective gravitational mass (m) of O [S]. The regularity 
conditions on x a-re 7] 



|x(0)|<oo, x(0) = 0. 



(5) 



Throughout we refer to the boundary (where the 
isotropic pressure vanishes and the solution matches onto 
vacuum) as E. Because of the form of ((T]), it is important 
to recall the following scaling property: If a spherically 
symmetric perfect fluid is represented by ds^ , then under 
the conformal transformation ds^ — > 5'^ds^, where (5 is a 
constant, p ^ pjS^, p — > p/S^ and m 6m. 



III. TOLMAN IV 



As an introduction to what follows, we start with a 
consideration of the Tolman IV solution @| • This so- 
lution is generated here in the following way: We choose 



1 



2„2 ■ 



(6) 



Since C2 can be absorbed into the scaling property men- 
tioned above and C3 can be absorbed into the scale of r, 
there is no loss in generality in setting C2 = C3 = 1. It 
now follows from ^ that 



M = 



'(l-Ci(l-r^)) 
l + r2 



(7) 



The choice ([6|) obviously introduces a coordinate singu- 
larity at r = 1. This is discussed below. The isotropic 
pressure and energy density are given by 



iirpir) — 



r2(l + 4Ci) + 2(l + Ci) 
1 + r2 



(8) 



2 



and 



87rp(r) 



r4(l+4Ci)+r2(5 + 2Ci) + 6Ci 



(l + r2)2 



(9) 



respectively. Bearing in mind the scale factor (5, we note 
that Ci = -47rp(0)/3 = 47rp(0) - 1. The adiabatic sound 
speed is independent of Ci and is given by 



1 + 

5-3r2 



Solving for Ci at S we obtain 



Ci 



2 + H 



2(l + 2r2) 



(10) 



(11) 



which shows that —1 < Ci < —1/2 and it follows that 
both p and p are monotone decreasing. Some bound- 
ary properties of this solution are shown in FIG. [1] and 
some internal properties are shown in FIG. [H The case 
Ci — —1/2 does not correspond to a finite object and 
the associated total mass m diverges as r — » 1 = E. 
This case is also distinguished by the peculiarity that 
all polynomial invariants of the Riemann tensor van- 
ish identically as r — > 1. The coordinate singularity at 
r — I plays no role throughout the range of finite objects, 
-1< Ci < -1/2. 
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FIG. 1: Some boundary properties of the Tolman IV solution. 
The abscissa is r^. defined by p{rj,) — 0. For the top curve 
the ordinate is —Ci. The next gives p and so gives the den- 
sity discontinuity at the boundary. The last curve gives the 
effective gravitational mass m, the total mass of each object. 
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FIG. 2: Some internal properties of the Tolman IV solution. 
The horizontal curve gives the adiabatic sound speed for all 
Ci. The thick curves give p/p{0) and the other curves give 
p/p(0). The values of Ci, right to left for p/p{0) and bottom 
up for p/p(0), are -0.5, -0.75, -0.85 and -0.95. 



the solution appears in the more familiar form 
(1-F2r2)dr2 



(l + (l + 2Ci)r2)(l-f r2) 



-I- - (1 + r^)dt\ 



(13) 



IV. INTEGER SOLUTIONS 



We now consider 



r2 

„2x(r) _ _ ^2 



(1-C|r2)^ 



(14) 



where n is and integer > 1. As explained above, there 
is no loss in generality in setting C2 = C3 = 1. It now 
follows from ^ that 

(_l)"-2^3 

^^= (l_;2)n-2p i^^rH + C) (15) 



where 



8 3n — 2 



Under the coordinate transformation 
Vl-r(r)2 



(12) 



{2n ~ l)r^ + 1, 



(17) 
(18) 



3 



and 



H = 2Fi{3-n, 



3n — 2 5n — 3 e 



2n - 1 ' 2n - 1 ' 2n 



(19) 



Here F is the standard hypergeometric function defined 

by 



2Fi{a,b;c;z) = ^ 



fe=0 



{a)k{b)k 
{c)k k\ 



(20) 



where {x)k is the Pochhammer symbol 

_ r(.T + k) 



T{x) 



= x{x + l)...{x + k + l). (21) 



All these solutions are physically distinct from, but bear 
some similarities to, the Tolman IV case ri = 1 9]. The 
case n = 2 is treated separately in the Appendix. For 
n > 2 the effective gravitational mass is given by 

M{2{n - Ifr'^ + 4(n - ly"^ + 2) - n{n - 2)r^ - 2m 



2(1 -r2)2+"/2 



where M is given by ((T3]). 
Let us define the "tenuity" 



A 



(22) 



(23) 



where A signifies the area of the two-surface t — const 
and r = r\^. Equivalently 



a = 



We find 



9ee 

(3?i - l)r2 + 1 



(24) 



(25) 



for all n and we note that a > 3 for all finite objects in 
this class of solutions. 



A. Limiting solutions 

As in the case of Tolman IV, for each ti > 1 there 
exists a global solution, characterized by a unique value 
of Ci. These do not correspond to a finite object and 
the total mass m diverges as r ^ 1 = E (as does the 
coordinate mass M for n> 2). These solutions share the 
distinction that all polynomial invariants of the Riemann 
tensor vanish identically as r — > 1. For n — 2, the critical 
value of Ci is 2^/^. For n > 2 the critical values are given 
by 



Ci = -ii^ Pr=inr=l) 



(26) 



where -0 is given by Vr=i and Hr=i by (fT7|) and 

respectively with e = 2n. Some properties of these 
global solutions are shown in FIG. [31 
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FIG. 3: Some internal properties of the global solutions. For 
the thin curves n — 2, for the thick curves n — 5 and for 
the dashdot curves n — 8. The horizontal curves give the 
adiabatic sound speed. The vertical curves give p/p{0) on top 
and p/ p(0) below. 



B. Even Integers 

For n even, all allowed values of Ci are > and for 
n > 2 the maximum allowed value of Ci is obtained from 
the limit r = T, and the minimum allowed value of 
Ci is obtained from the limit r ^ 1 = E as given by (|26p . 
For n = 2, 2^/^ < Ci < 4 and for n>2 

- (V P,.=iHr=i) < Ci < -(0 Hr^o) + 2n (27) 



where Tir=o is given by (|T9)) with e = 1. The coordinate 
singularity at r = 1 plays no role throughout the range 
of finite objects. Bearing in mind the scale factor 6, we 
find 



Ci = -47rp(0) - (V' Hr^o) + 2n 
and note that 



p(0) 



P(0) 



47r 



(28) 



(29) 



as in the case of Tolman IV {n — 1). Some boundary 
properties of the solutions for even n are shown in FIG. 2] 
and some internal properties are shown in FIG. [5l 



C. Odd Integers 

For n odd, all allowed values of Ci are < and the 
minimum allowed value of Ci is obtained from the limit 
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FIG. 4: Some boundary properties for even integers. For 
the thin curves n — 2, for the thick curves n — 4 and for 
the dashdot curves n = 6. The abscissa is defined by 
p(rj,) = 0. For the top curve the ordinate is Ci scaled by its 
maximum value at the center. The next gives p/p(0) and so 
gives the scaled density discontinuity at the boundary. The 
last curve gives the effective gravitational mass m, the total 
mass of each object. 




FIG. 5: Some internal properties of the n = 2 solution 
(explicit properties given in the Appendix). The horizon- 
tal curves give the adiabatic sound speed. The thick curves 
give p/p{0) and the other curves give p/p(0). The values of 
Ci, right to left for p/p{0) and bottom up for p/p{0) and the 
sound speed, are 2^'''^, 3.5, 3.75 and 3.9. Solutions with n > 2 
show qualitatively the same behaviour. 



r — > = S and the maximum allowed value of Ci is 
obtained from the limit r ^ 1 = S as given by (|26p . For 
n odd we have 



(30) 



Again, the coordinate singularity at r = 1 plays no role 
throughout the range of finite objects. Again, bearing in 
mind the scale factor S, we find 



Cl = 47rp(0) - {ip Hr^o) ~ 2n 



(31) 



and also holds. Boundary properties of the solutions 
for odd 71 are qualitatively similar to FIG. [4] and their 
internal properties are qualitatively similar to FIG. [51 



is the Tolman IV case n — 1. In order to transform the 
foregoing solutions into explicit isotropic coordinates, we 
must solve for r(r) from 



ln{r) 



dr 



ry/1 - 2M{r)/r 



(33) 



where M{r) is given by (jlSp. We have found that this is 
not possible for any n. 



NON INTEGER SOLUTIONS 



Let us now consider (|14|) but with n generated by 



D. Coordinate Transformations 

In order to transform the foregoing solutions into ex- 
plicit curvature coordinates, we must solve for r(r) from 
the generalized form of (fT^ . 



(32) 



^(r) 



(1 - r(r)2)"/2 



Whereas this is possible in principle only for n < 5, the 
only value of n which we find gives a manageable metric 



l + 2i 



(34) 



where i is an integer > 1. This choice also admits global 
solutions but only with Ci = 0. These solutions all have 
monotone decreasing p but p fails to be monotone (in- 
deed p(0) = for z = 1) and so these solutions have to 
be considered unphysical. Some properties of these solu- 
tions are shown in FIG. [S]for completeness. Non global 
solutions exist for Ci 7^ but for these p also fails to be 
monotone. 



5 



1 . 6- 



1 . 4- 



recherche sur la nature et les technologies. PJE is sup- 
ported by a scholarship from the National Science and 
Engineering Research Council of Canada. KL is sup- 
ported by a grant from the Natural Sciences and Engi- 
neering Research Council of Canada. Portions of this 
work were made possible by use of GRTensorll p^ . 



1 . 2- 



APPENDIX: n = 2 
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For n = 2 we find 
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M 



(l + 3r2)4/3' 



(A.l) 



. 2- 



. 2 
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r^2{l + 3r'^)^/^ -Cijl+r'^f) 
(l-r2)3(l + 3r2)4/3 ■ 



(A.2) 



FIG. 6: Some properties of the non integer global solutions. 
For i = 1 (dashdot curves) p andp/p(0) are shown. The other 
cases show p/p(0) and p/p{0). In all cases p/p(0) is monotone 
decreasing but p/p(0) is not. The values of i are 2, 3, 4, 5 
and 100 (essentially the same as n = 1/2). These are read 
top down on both sets of curves. 



2(1 + 3r2)4/3(2 + r2) - Ci(l + 5r2)(l + r^) 

^ ' (l + 3r2)4/3 ' 

(A.3) 

so that 



VI. DISCUSSION Ci ^ -47rp{Q) + A (A.4) 

The Einstein static Universe has been transformed into and finally 
an explicit and exact two - fold infinity of physically ac- 
ceptable perfect fluid solutions of Einstein's equations. 6 4 2 27/3 
The solutions are physically interesting because they are _ 2 ^^^^^ + 29r +117' -|-3) — 6(l + 3r ) 
useful for the study of internal properties of neutron (1 + 3r2)''/3 
stars. The emphasis here is on the role of coordinates. 
Whereas the entire class of solutions can be considered so that 
as generalizations of the familiar Tolman IV solution, no 

member of the class can be written explicitly in isotropic 4 

coordinates and, except for a set of measure zero, no ^1 ~ '^^Pi^) + 2. (^-6) 

member of the class can be written explicitly in curva- 
ture coordinates either. „, rii 1- 1 i- 1 1 • • 1 

ihe square ot the adiabatic sound speed is given by 
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